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ABSTRACT 

We have measured the time delay between the two images of the gravitationally lensed quasar HE1 104—1805 
by combining observations made with the Wise Observatory lm telescope and published observations of this 
system by Schechter et al., spanning a total of five years, from 1997 to 2002. Based on a variety of techniques, 
we find that the best fit time delay is — 161^Z'^jj days. The 68% and 95% confidence intervals include the 
systematic errors due to an observed component of uncorrelated variability between images. The delay is 
shorter than predicted by simple models and may indicate a large external shear or a large value of the Hubble 
parameter, h > 0.75 (95% confidence). The residual light curve between the two images shows a longterm 
trend of ~ 0.04 mag yr _1 , presumably the result of microlensing by stars in the lens galaxy, but also short 
timescale (~ 1 month) variability with a mean amplitude of about 0.07 mag. 

Subject headings: cosmology: gravitational lensing — Quasars: general — Quasars: individual 
(HE1 104-1805) 
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1. INTRODUCTION 

Time-delay measurements in a large sample of gravitation- 
ally lensed systems could be used to measure the Hubble pa- 
rameter, Hq, directly at high redshift, therefore avoiding the 
systematics due to distance ladder calibrations and large-scale 
motions (e.g., Turner, Cen, & Ostriker 1992). Alternatively, 
using a measured value of the Hubble parameter from other 
techniques (e.g., Bennett et al. 2003), measurement of a time 
delay allows discriminating between different mass profiles 
of galaxies (e.g., Kochanek 2002). In some systems, complex 
geometries caused by multiple lensing galaxies hinder the use 
of time delay measurements for Ho determination or for mass 
profile studies. Even in such systems, time delay measure- 
ments are important for isolating the microlensing variability 
from the intrinsic variability. Microlensing in lensed quasars 
could allow an independent measurement of the fraction of 
dark matter in galaxies (Schechter & Wambsganss 2002), and 
can be used to put limits on the size of the continuum emission 
regions in quasars (e.g., Wambsganss, Paczynski, & Schnei- 
der 1990; Wambsganss 2002; Wyithe, Agol, & Fluke 2002). 

HE1 104— 1805 is a double-image lensed quasar discov- 
ered by Wisotzki et al. (1993). The image separation is 
A0 = 3.19", the source redshift is z, = 2.319 (Wisotzki et 
al. 1993), and the lens redshift is zi — 0.729 (Lidman et al. 
2000). This is an unusual system, in which the image closer 
to the lens (image A) is the brighter one. Wisotzki et al. (1993) 
noted the emission line flux ratio between the images is 2.8, 
and that component A has a bluer continuum and lower equiv- 
alent width emission lines than B. They interpreted the dif- 
ferences as being caused by microlensing of the quasar con- 
tinuum source in the brighter image. Wisotzki et al. (1995) 
reported that the continuum flux in both images is highly vari- 
able, but that the line fluxes do not change. Again, this could 
be interpreted as evidence for strong microlensing. Lehai et 
al. (2000) have attempted modeling this system. They as- 
sumed an image flux ratio of about 4, and found this limits 
the models to a narrow range of predicted time delays. 

Wisotzki et al. (1998) measured the continuum flux in 19 
low resolution spectra, taken with the ESO 3.6m telescope 
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between 1993 to 1998. By matching the light curves of the 
two images, they favored a time delay of about —270 days 
(in the sense that the B image leads the A image), although 
they concluded that a value as short as —100 days could not 
be excluded. Gil-Merino, Wisotzki, & Wambsganss (2002) 
analyzed the B-band acquisition frames for the spectroscopy 
of Wisotzki et al. (1998) and found a best fit time delay of 
-310 ±20 days (2(7 errors). Pelt, Refsdal, & Stabell (2002) 
re-analyzed the Gil-Merino et al. (2002) observations, and 
argued that the time delay is somewhere between —330 and 
—255 days, but with a large uncertainty. In a recent paper, 
Schechter et al. (2003) presented 3 years of V-band photome- 
try of HE1 104-1805, obtained on 102 nights with the OGLE 
1.3m telescope. They did not succeed in finding a consis- 
tent time delay for the system. In fact, the root-mean-square 
(rms) of the difference between the component light curves 
as a function of time delay showed two minima, at about 
— 150 days and —360 days. Moreover, Schechter et al. (2003) 
noted that the structure functions of images A and B are quite 
different, with image A being more than twice as variable as 
image B on timescales of less than a month. Thus, there is 
presently no clear evidence of a time delay in this system, let 
alone an unambiguous measurement of its value. 

In this paper we present new photometric data for 
HE1 104— 1805 from the Wise Observatory lens monitor- 
ing project. Combined with the OGLE observations of 
HE1 104-1805 from Schechter et al. (2003), our data es- 
tablish the delay securely and isolate the various correlated 
and uncorrelated variability components. We will assume 
throughout a cosmology with £l m = 0.3, CIa — 0.7, and Hq = 
100/; km s- 1 Mpc- 1 . 

2. OBSERVATIONS AND REDUCTION 

For the past three years we have been monitoring a sam- 
ple of 27 known gravitationally lensed quasars and lens can- 
didates with the Wise Observatory lm telescope. The ob- 
jects are monitored on a weekly basis in the Johnson-Cousins 
/?-band, with occasional coverage in the /-band or V-band. 
Frames are obtained with a cryogenically cooled Tektronix 
1024 x 1024-pixel back-illuminated CCD. The scale is 0.7" 
per pixel. The median seeing of about 2" does not al- 
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low resolving most of the lensed objects. However, some 
of the pairs are resolved (e.g., Q0142-100, HE1104-1805, 
RXJ0921+4528, Q0957+561). The data are reduced auto- 
matically on a daily basis. In each field we perform aperture 
photometry of all sources in the frame and the preliminary 
combined-image light curves are posted on the WWW 1 . 

We present here 7?-band measurements of HE1 104— 1805 
taken between 1999, November 14 and 2002, June 18. The 
combined OGLE and Wise observations span over five years, 
from 1997 to 2002. The Wise observatory photometry of 
HE1 104— 1805 was performed using a special program we 
have written for this purpose. All the frames are aligned and 
rotated to the same reference frame. The program constructs 
a numerical point-spread function (PSF) from preselected ref- 
erence stars. Given accurate positions of the lensed images 
from Hubble Space Telescope images (Lehar et al. 2000), 
it fits the PSF to the images and minimizes the % 2 (assum- 
ing Poisson noise) with respect to the two free parameters, 
the magnitude of the two images, m\ and mi. In the case of 
HE1 104-1805, the PSF was constructed from 12 stars. After 
rejecting all measurements with errors larger than 0.1 mag, 
we are left with 79 epochs for image A and 49 epochs for im- 
age B 2 . Note that the light contributed by the lensing galaxy 
to image A was neglected in the photometry, as it contributes 
only 0.3% and 3.5% of the total light of image A in the V and 
/ bands, respectively. 

The error bars were calculated from the covariance matrix 
of the x 2 surface (e.g., Press et al. 1992). To validate the relia- 
bility of the error bars, we have made an empirical, magnitude 
dependent, error estimate for each lensed image and com- 
pared it with the x 2 error estimate. The empirical errors were 
calculated as follows. Using the same algorithm, we have 
measured the difference, Am', for each frame taken at epoch 
t , between the magnitude of each reference star, i, and its av- 
erage magnitude over time, m* (with the magnitudes measured 
relative to the 12 co-added reference stars). The error in the 
measurement of a particular source will depend on its intrinsic 
magnitude and on the observing conditions on the particular 
epoch. To find this dependence, we fit a function, with two 
free parameters p and q, of the form: log Am\ = p + qihi to all 
the frames, and find q = 0.3. This function approximately de- 
scribes the expected dependence of error on magnitude. Note 
that, in the case that the background Poisson noise is the dom- 
inant source of error, we expect q — 0.4, while in the limit of 
zero background we expect q = 0.2. Assuming that q = 0.3 is 
constant over all nights, for each individual frame we refit this 
formula to obtain p t . This fit gives, for each frame, the typical 
error as a function of magnitude. Finally, for each frame, we 
have compared the x 2 errors in the magnitudes of the lensed 

images, A and B (mf, mf ), with the empirical errors \QP< + ^ 
and IQP'+i" 1 ? . We find excellent agreement between the two 
error estimators. 

Table 1 lists the reference star coordinates, magnitudes, m,, 
and the 68% confidence interval of their magnitude distribu- 
tion (68% CI). All magnitudes in this paper are given relative 
to reference star ST. 

3. WISE-OGLE INTERCALIBRATION 

We have intercalibrated the OGLE and Wise photometry, 
in order to produce combined light curves with a long time 

1 http://wise-obs.tau.ac.il/~eran/LM/ 

2 The measurements, as well as finding chart with reference stars, are avail- 
able from http://wise-obs.tau.ac.il/~eran/LM/HE1104/ 



Table 1. PSF reference stars 



Ref R.A. (J2000) Dec. m, 68%CI 



SI 


11:06:33.17 


-18:21:39.1 


0.262 


0.025 


S2 


11:06:33.86 


-18:20:29.6 


0.585 


0.038 


S3 


11:06:44.75 


-18:20:43.2 


0.374 


0.035 


S4 


11:06:38.06 


-18:19:10.7 


0.344 


0.043 


S5 


11:06:30.38 


-18:19:33.0 


1.633 


0.088 


S6 


11:06:29.74 


-18:20:02.3 


1.321 


0.065 


SI 


11:06:23.23 


-18:19:41.5 


0.000 


0.028 


ss 


11:06:22.50 


-18:21:30.4 


0.384 


0.038 


S9 


11:06:35.63 


-18:21:29.3 


0.987 


0.055 


SW 


11:06:37.00 


-18:22:45.5 


0.414 


0.033 


Sll 


11:06:37.39 


-18:23:00.3 


1.551 


0.100 


512 


11:06:38.31 


-18:22:48.1 


1.284 


0.060 



NOTE. — Astrometry is based on 38 USNO-A2.0 stars (Monet 1998) with 
0."4 rms in each axis. Magnitudes are given relative to reference star 57. 

coverage. To this end, we take advantage of the fact that 
the Wise and OGLE light curves overlap in the year 2000 
(J.D. - 2451450-2451750). However, the OGLE observa- 
tions were obtained in the V band, while the Wise observa- 
tions were obtained in the R band. Wisotzki et al. (1995) have 
already noted that the flux ratio between the quasar images 
depends on wavelength. We therefore intercalibrated the pho- 
tometry for the light curve of each lensed image separately. 
For each image, we performed a x 2 fit between the overlap- 
ping OGLE and Wise light curves to determine the magnitude 
offset and an optional stretch factor. We allow for a stretch 
factor since quasars are known to vary in color, and the ampli- 
tude of variation in the different bands is thus not necessarily 
the same. We used two types of fits: (i) An offset + stretch 
factor, niR = a + bmy, where m,R and my are the magnitudes 
in the R and V bands, respectively; and (ii) A pure magnitude 
offset, niR = a + my- 

Since the OGLE and Wise observations during the overlap 
period were not carried out on the same days, we linearly in- 
terpolated the OGLE observations (which are more frequent 
and have smaller error bars) to the times of the Wise obser- 
vations. We set the error of each interpolated magnitude to 
be 

£ = ^£ 2 + sf{At min ), (1) 

where £„ is the error of the real measurement nearest in time 
to the interpolated point, s/(Af, m „) is the value of the struc- 
ture function of the interpolated light curve at lag Af min , and 
At m i n is the time between the interpolated epoch and the near- 
est epoch which has a real measurement. The structure func- 
tion is defined as 

sf 2 (t)=2[V-DCF(t)}, (2) 

where V is the variance of the light curve and DCF{t ) is the 
(unnormalized) discrete auto-correlation function (Edelson & 
Krolik 1988), calculated with 14-day bins. This scheme gives 
the interpolated points a realistic weight in the x 2 minimiza- 
tion. 

The 1 a errors of the pure magnitude-shift fit between the V 
and R light curves are less than 0.03 mag for both images. Ap- 
plying the shift+stretch fit, the best-fit stretch factor between 
the V and R light curves is consistent with unity and does not 
improve the fit significantly. As we will show below, allowing 
for a non-unity stretch factor does not change the best fit time 
delay. Therefore, in what follows we will use the magnitude- 
shifted light curves (i.e., with no stretch). We test in |]4]for the 
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1 $50 1500 1550 1600 1650 1700 1750 

JD - 2450000 

FIG. 1. — Image A (upper panel) and image B (lower panel) light curves 
during the period of overlap between the V-band OGLE observations (tri- 
angles) and the R-band Wise observations (circles). The observations were 
brought to the same scale by % 2 minimization. Magnitudes are relative to 
reference star 57. 



Table 2. Intercalibrated image A light curve 



JD-2450000 


mag 


error Observatory 11 


666.46490 


-0.418 


0.005 1 


666.47390 


-0.418 


0.005 1 


760.84802 


-0.360 


0.003 1 


760.85665 


-0.354 


0.003 1 


762.85473 


-0.341 


0.004 1 



NOTE. — Table 2a in its entirety, as well as Table 2b for image B, are 
available via the electronic version. 

a l- Observations from Schechter et al. (2003; OGLE V band); 2- Our ob- 
servations (Wise R band). 



uncertainty in the time delay induced by the range of possible 
shift+stretch parameters. 

The combined light curves in the zone of overlap are shown 
in Figure [2 for images A (upper panel) and B (lower panel). 
In Figure [2 we show the complete (1997 -2002) light curves of 
images A (upper panel) and B (lower panel), after application 
of the intercalibration to the whole dataset. The complete 
intercalibrated data, from Schechter et al. (2003) and our ob- 
servations, are given in Table 2. 

4. TIME DELAY BETWEEN IMAGES A AND B 

We now search for the time delay between images A and 
B by using x 2 minimization. From Figure |2] it is apparent 
that the long term variation of image A is more pronounced 
than that of image B. This could be the result of microlens- 
ing in the system. Thus, in the % 2 minimization we leave as 
free parameters not only a time delay, but also a linear trend 
of magnitude with time between the light curves. The fit is 



described by 

Y 2 = y [mf + ,-mf+S{t-t mid )-C] 2 

X V e%+sf(At min ) ( ' 

where mf and mf are the magnitudes as a function of the time, 
f, for images A and B, respectively. The arbitrary constant f m ,y 
is defined as the midpoint between the first and last observa- 
tions (i.e., JD — 2451525.378), C is the mean magnitude dif- 
ference between the images, and S and z are the fit parameters 

— a linear trend between the light curves, and a time delay, re- 
spectively. The fit was performed, simultaneously for S and C, 
between the light curve of image A and each of the T-shifted 
light curves of image B. In this process, the light curve of 
image A at epoch t was interpolated to the times (f + t) of 
the image B observations. Again, we used linear interpolation 
and set the error of the interpolated magnitude using Eq. \l\ 
The bin size used for calculating the DCF necessary for de- 
termining the structure function, sf (At „„■„), has a small effect 
(< 2 days) on the best-fit time delay. 

Figure [3] shows the % 2 P er degree of freedom (dof) of 
the light curve fitting as a function of the time delay. The 
dashed curve shows the dof. For each time delay, the best-fit 
linear trend (S) and magnitude offset (C) were recalculated. 
There is a distinct minimum in the % 2 at a time delay of 

— 161 day. The best fit parameters are S = 0.043 mag yr _1 
and t = -160.9 day (with x 2 /dof = 303/229 = 1.32). For 
these parameters, the magnitude difference (C) between the 
images in the R band is Atur = 1 .595 ± 0.004 mag. Figure|4] 
shows the light curve of image A (filled circles), overlayed 
by the slope corrected (S = 0.043 mag yr _1 ) and time-delay 
shifted (t = —160.9 day) light curve of image B (empty cir- 
cles). The formal % 2 errors on the time delay are less than 
1%. However, these errors are not realistic. As we will show 
below, there are several sources of systematic errors that need 
to be taken into account. 

As noted above, a systematic error is introduced to the time 
delay by the uncertainty in the OGLE and Wise intercali- 
brations (using the shift algorithm or the shift+stretch algo- 
rithm). We have used the uncertainty in the shift+stretch pa- 
rameters and their covariance term, and have run 1000 Monte- 
Carlo simulations in which we randomly drew shift+stretch 
parameters from a bivariate Gaussian distribution, and re- 
peated the time-delay fit. We find a time delay distribution of 

— 159.2^j'^2 3 n day (68% and 95% confidence errors). The 
short time delays, found in some of these simulations, of 
about — 1 30 day, are the result of stretch-factor values smaller 
than b = 0.75. Quasars are known to have larger variability 
amplitudes at bluer wavelengths, corresponding to b < 1 . For 
example, Giveon et al. (1999) found A(B - R) w 0.25AB for 
a sample of quasars at z ~ 0.2, which corresponds to a stretch 
factor b = 0.75 between B and R. For HE1 104-1805, we are 
observing a smaller wavelength interval, between V and R, but 
on the other hand, a higher redshift. The V and R bands sam- 
ple the restframe UV, which may have a stronger dependence 
of variation amplitude on wavelength. Given the current un- 
certainty in the stretch factor, there is about 4% probability 
that the time delay between the images of HE1 104— 1805 is 
between — 125 and — 155 days. If we limit the stretch parame- 
ter to 0.75 < b < 1, then the 95% confidence error on the time 
delay, due to the stretch + shift uncertainty, is smaller than 
2 days. 

With these results we can reject the possibility that the time 
delay is in the region of —250 to —330 days, as suggested 
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by Gil-Merino et al. (2002) and Pelt et al. (2002). Given 
the significant microlensing variability observed in this sys- 
tem (e.g., Schechter et al. 2003), measuring the time delay 
from sparsely sampled light curve is difficult. 

Figure |4] shows that the main features and trends appear in 
the light curves of both images. However, there are significant 
coherent fluctuations in the difference light curve, as already 
noted by Schechter et al. (2003). Figure [5] shows the dif- 
ference between the A and B light curves after applying the 
best-fit linear trend and time delay. In order to subtract the 
two light curves we interpolated the time-delay- and linear- 
trend-corrected light curve of image B to the times of image- 
A observations. Note that this plot mixes the V- and /?-band 
measurements. The error bars were calculated according to: 

Am = ^ej + el n + sf(At mi „), (4) 

where Ea is the error for image A, and Eb„ is the error for 
image B at the point nearest to the interpolation time. The 



structure function of the residual light curve, shown in Fig- 
ure[6] rises rapidly from to 50 days lag (time in the observer 
system) and then stays approximately constant at a level of 
5 x 10~ 3 mag 2 (w 0.07 mag). The structure function was 
calculated using Eq. |2j where the DCF was calculated with 
14 days bins. 

The short-timescale, uncorrelated, variability in the light 
curves is a second source of systematic error in the delay mea- 
surement. To evaluate the effect of these fluctuations on the 
delay accuracy, and to test the hypothesis that the time delay 
might be completely different (e.g., Gil-Merino et al. 2002), 
we have performed additional Monte-Carlo simulations. We 
have used the algorithm of Timmer & Koenig (1995) to ar- 
tificially generate random red-noise light curves. The light 
curves were generated with the same temporal sampling pat- 
tern as the real observations obtained by OGLE and by us, and 
with a power density spectrum (PDS) proportional to f~ a , 
where / is the frequency and a is the power-law index. The 
PDS of quasar variability in the optical range is poorly known 
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-400 -300 -200 -100 100 200 300 400 6 ° 

Time delay (days) 

FIG. 3. — The % 2 /dof as a function of delay time (solid curve). For each 
time delay, the best-fit linear trend (S) and magnitude offset (C) were used. 
The dashed curve shows the dof. 



(see Markowitz et al. 2002, for the PDS in X-rays) but is 
likely a power law with a between 1 and 2 (Giveon et al. 
1999). The standard deviation (StD) of the variations about 
the mean of the simulated light curve was scaled to the StD 
of the observed image-B light curve. To simulate the image- 
A light curve, we duplicated the image-B light curve with a 
time delay Tj. We then added to each point in both simulated 
light curves a normally distributed noise with StD taken from 
the individual errors in the observed light curves. In order to 
mimic the uncorrelated noise, we also added, in some of the 
simulations, the residual light curve (shown in Fig. |5j to the 
simulated image-A light curve. We then used % 2 minimiza- 
tion, as for the real data, to search for the best time delay. For 
each set of parameters (e.g., a, %s) we repeated this process 
300 times. Note that we do not take into account the linear 
trend in the simulations. 

Our main conclusions from the simulations are: (i) The 
power-law index a has the largest effect on the uncertainty 
of the deduced time delay, with a ~ ±1 day 95% -confidence 
error for a = 2, and a ~ ±7 day 95% -confidence error for 
a = 1 . A small value of a introduces variability on timescales 
shorter than the mean sampling interval; (ii) When the uncor- 
related variability is added to the simulated light curve of im- 
age A, the confidence interval increases by about 30%; (iii) 
Testing for T s 7^ —161 day (e.g., t s « —300 day; as suggested 
by Gil-Merino et al. 2002) gives similar uncertainties in the 
measured time delay to those found using t s = —161 day. 
Combining the uncertainty from these simulations (a = 1), 
and the systematic errors described above, we have adopted a 
time delay of — day (68% and 95% -confidence er- 

rors). Moreover, we can reject with high confidence the pos- 
sibility that the time delay is in the vicinity of ~ —300 day, as 
suggested by earlier works on HE1 104—1805. 

To verify our result using another method, we have applied 
the ZDCF cross-correlation technique (Alexander 1997) di- 
rectly to the data. We find a best fit time delay of about 
— 156^ 5 day, with a peak correlation of 0.82 ±0.03, con- 
sistent with our X 2 fit- The result changes by less than 2 days 
if, before cross-correlating, we subtract the linear trend (S) 
from the image-A light curve, subtract first degree polynomi- 



als from both light curve, or use the image fluxes instead of 
magnitudes. 

5. DISCUSSION AND SUMMARY 

We now discuss briefly the implications of the time delay 
we have found, T = —161 ±7 days. Lehar et al. (2000) mod- 
eled the lensing potential of HE1 104— 1805 using a singular 
isothermal ellipsoid or a constant mass to light ratio, with and 
without the external shear expected from galaxies projected 
within ~ 20", assuming they are at the lens redshift. Their 
predicted time delays for these four different models are in 
the range of — 129/i 1 to —263h~ l days. Assuming h = 0.7 
(Bennett et al. 2003), the shortest time delay (—184 days) 
predicted by the singular isothermal ellipse + external shear 
model is larger than our result. Note that HE1 104— 1805 is 
a double-image lensed quasar, and therefore it was necessary 
to use the flux ratio between images in order to constrain the 
lens model. Interestingly, Lehar et al. (2000) note that an 
external shear twice as large as their most extreme model is 
needed in order for the lens mass to be aligned with its light. A 
larger shear would lower the predicted time delay, and make it 
more consistent with our measurement and with current mea- 
surement of Hq. Alternatively, if we adopt the model range of 
Lehar et al. (2000), our measured limit of T > — 172 day (95% 
CL) sets a lower limit of h > 0.75. HE1 104-1805 apparently 
is an unusual system, as indicated by the fact that the bright 
image A is the one closer to the lens and that the measured 
time delay is not consistent with the time delay predicted by 
simple models. We defer dealing in more detail with the im- 
plications of the time delay to a future work. 

The mean flux ratio of A /B ~ 4.4 we find in this work 
is significantly different from the emission line flux ratio of 
A/B ~ 2.8 reported by Wisotzki et al. (1993). Extrapo- 
lating the linear trend we have found, S ~ 0.04 mag yr , 
into the future, suggests that the broad-band flux ratio 
between the images will decrease to the level of the 
emission-line flux ratio (A/B ~ 2.8) in about a decade. 
At the source and lens redshifts of HE 1104-1805, the 
Einstein-radius crossing time for stellar objects in the 
lens galaxy, having mass M and transverse velocity v is 
2O.3(M/M0) 1 / 2 (v/6OOkms" 1 )- 1 (/i/O.7)- 1 /2 yr . Thus, the 
slow trend is well explained by microlensing in the macro im- 
ages of HE1 104- 1805. 

As shown in Figures |4] and |5] even after removing a 
longterm linear trend, there is significant uncorrelated vari- 
ability. As already noted by Schechter et al. (2003), judging 
from the amplitude and timescale of the variability, it seems 
that most of the uncorrelated variability occurs in image A (the 
one nearest to the lens). Uncorrelated short-timescale vari- 
ability has been observed in other lensed system (e.g., Burud 
et al. 2002). However, in HE1 104-1805 it has larger ampli- 
tudes, with a mean of 0.07 mag (up to 30% peak-to-peak) on 
timescales of less than a month. 

Wambsganss et al. (1990) have pointed out that, for large 
microlensing optical depth, some events can have timescales 
that are considerably shorter than the Einstein-radius cross- 
ing time. This is caused, for example, by the passage of the 
source near a cusp. Wyithe & Loeb (2002) have suggested 
that the low-amplitude, fast (~ 50 days), uncorrelated vari- 
ability observed in RXJ0911+05 (Hjorth et al. 2002) and 
in SBS 1520+530 (Burud et al. 2002) can be explained by 
stellar microlensing of a smooth accretion disk that is oc- 
culted by optically thick broad-line clouds. Although their 
model predicts the short timescale variability, it cannot repro- 
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FIG. 5. — The difference light curve between images A andB, after applying 
the best-fit linear trend and time delay. 
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FIG. 6. — Structure function of the difference light curve of Fig.l5l 



duce the large amplitude observed in HE1 104— 1805. An- 
other possibility they consider is microlensing by planetary- 
sized objects of 10~ 2 — 10~ 4 M , but they show that this 
cannot produce variations with 10% amplitude or greater. 
Gould & Miralda-Escude (1997) have suggested that the mi- 
crolensing of hot spots (or any structure) in the fastly rotat- 



ing quasar accretion disk can give rise to fast variations with 
timescales of ~ 1 month. Wyithe & Loeb (2002) have simu- 
lated microlensing of an accretion disk with 100 spots. They 
find a typical variability amplitude of ~ 10% on one-month 
timescales. Schechter et al. (2003) have argued that the vari- 
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ability in HE1 104— 1805 is best explained by such a model 
with v/c ~ 0.25. In order to test these models, we currently 
continue to observe this system frequently in the VRI bands. 
With more data at hand, we will attempt to re-address the na- 
ture of the microlensing in this system in a future paper. 

To summarize our main results, we have measured the time 
delay of the lensed double-image quasar HE1 104— 1805. We 
have used detailed simulations at every stage of the reduc- 
tion and analysis in order to obtain realistic error estimates. 
We find that the light curves are best fit with a time delay of 
— 161 ±7 days, and with a linear trend between the images of 
about ~ 0.04 mag yr . Our measurements resolve the pre- 
vious ambiguities pertaining to the time delay in this system. 
However, the time delay is shorter than predicted by any ex- 



isting models. The linear trend is likely due to stellar-mass 
microlensing in the lens galaxy. We confirm previous reports 
that the residual light curve between the linear-trend-corrected 
and time-delay-shifted light curves shows significant variabil- 
ity on short timescales of about one month. Multi-band and 
spectroscopic observations could help explain the nature of 
this fast variability. 
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